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SUMMARY 

A theoretical treatment of absorption problems is presented 
in which the following cases are simultaneously considered: 

(a) Radiation is normal to an absorber of which the stations 
are plane parallel surfaces. 

(b) Radiations are of several polyenergetic types. 

(c) Induced radioactive isotopes decay to stable atoms in 
multistep decay processes. 

(d) Radiations from the absorber affect the time -dependency 
of the source activity. 


. INTRODUCTION 

Matrix methods have been applied to the problem of calculating 
transmission, reflection, and conversion to heat of time -dependent 
nuclear radiations (references 1 and 2). The analyses that were 
made, however, assumed single-step decay of induced radioactive 
atoms to stable ones. Many isotopes decay in several steps to 
stable materials, emitting a variety of types of radiation during 
the process. The trans uranic radioactive Isotopes are examples of 
multistep decay elements, as are certain isotopes of elements 
commonly used as structural materials such as iron and copper. 

The matrix methods of references 1 and 2 developed at the 
NACA Lewis laboratory were therefore extended to include the multi- 
step decay processes. The addition of these processes to the prob- 
lems considered In references 1 and 2 makes the defining equations 
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quite lengthy even though matrix notation is used. In order to 
facilitate presentation, tensor notation and the summation conven- 
tion, when feasible,- are used. 

The ass umpt ion is made in references 1 and 2 that the source 
intensity is unaffected by back-reflected radiation from the 
absorber; that is, the time -dependency of the source is assumed 
to be known explicitly for all times of interest. In practice, 
however, the radiation back-reflected to the source includes radi- 
ation from the radioactivities within the absorber. Thus, although 
it is unlikely that the source will be noticeably affected by radi- 
ation reflected in the literal sense, enough of the aforementioned 
radiation originating within the absorber may quite possibly reach 
the source to induce a new radioactivity therein. If the contri- 
bution of the new activity to the total source intensity is not 
negligible, the time-dependecy of the radiation from the source 
will be significantly altered. 

The change in source intensity as a function of time des- 
cribed in the preceding paragraph can be quantitatively calculated 
by methods shown herein. 


SYMBOLS AND NOTATION 


The following symbols are used In this report; 


[a] > [ B ] ) L C 3 > D3 ma ' trix coefficients 
E energy 


H 


thermal power generated 


h 

I 


rate of conversion of energy of nuclear radiation 
absorbed to thermal energy 

power of radiation, going to right, incident upon 
or emerging from station of absorber 


N number of radioactive atoms 

n number of stations in absorber 

one -half of power of radiation from radioactivity 
emerging from station of absorber 


P 
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q number of generations of daughters 

B power of radiation, going to left, incident upon 

or emerging from station of absorber 

r power. back-scattering coefficient 

t power -transmission coefficient 

K rate of con-version of radioactive energy of parents 

to radioactive energy of daughters 

A decay constant 

H rate of conversion of radiation absorbed to energy 

of radioactivity 

T time 

Superscripts and subscripts : 

( ( ) of energy e and type a in i^ station 

that becomes d tix generation radiation of 
energy e and type s 

8 8^k generation of daughters 

0 initial condition 

The following conventions have been adopted in the notation 

used; 


(1) If an index appears both as a subscript and as a super- 
script in the same term, a summation over the range of the index 
is indicated. As used herein, "tern" means any algebraic quantity 
or quantities separated from other algebraic quantities by +, -, 
or e signs. 

(2) In older to avoid confusion with the summation convention in 
tensor analysis, the further provision is made that only repeated 
Greek symbols will indicate summation. This t provision is made 
because English symbols may be repeated in the same tem but never 
both as a subscript and a superscript and will therefore not indi- 
cate summation. 
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(3) A symbol like q denotes the rate of conversion of 
radiation of energy e and type a absorbed in the i^ 11 station 
to radioactivity of a parent atom emitting radiation of energy e 
and- type s. Transformations from energy e and type a to 
energy e and type s will be symbolically indicated as in the 
preceding example. 


' ANALYSIS 

The analysis is based on the following assumptions : 

(1) Stations of the absorber are plane parallel surfaces . 

(2) Radiation is normal to absorber. 

(3) Part of the radiation absorbed at each station of the 
absorber is transformed to thermal energy and part to induced 
radioactivity in the absorber. 

(4) Radioactive atoms in the absorber decay in a finite num- 
ber of steps to stable atoms. 

(5) One-half of the radioactivity produced in each station of 
the absorber is emitted from each side of the absorber. 

(6) At seme instant of time, the time -dependency of the inci- 
dent radiation is known explicitly; 

Method . - Diagram (a) is included to clarify the discussion 
in the following paragraph. 


ii f 



ji+1 


pl 

pl 


R 1 



R i+1 







H 1 



As in references 1 and 2, 
the absorber is considered as 
being arbitrarily divided into 
n convenient stations normal 
to the path of the radiation. 
Radiation of power I 1 from 
the (i-l)^ station is inci- 
dent upon the i^ station of 
the absorber; whereas radia- 
tion R^- from the i^k sta- 
tion is incident upon the 
(i-l)'fck station. 
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In the treatment presented herein, the values of the following . 
physical constants most he known: 

h rates of conversion of energy of nuclear radiation absorbed to 
thermal energy 

r power hack-scattering coefficients 
t power -transmission coefficients 

K rates of conversion of radioactive energy of parents to radio- 
active energy of daughters 

A radioactive decay constants 

H rates of conversion of radiation absorbed to energy of radio- 
activity 

In addition, the following boundary conditions must be known 
for a complete solution to the problem* 

initial thermal power being emitted from i tJl station 
of absorber 

initial radiation power incident upon absorber from 
source 

one-half of initial power of radioactivity in 
i^ station of absorber for all i 

radiation power incident upon absorber from side of 
absorber opposite source 

The problem is to find I*, R 3 -, and H 3 - in terms of known 

quantities. The defining equations incorporating the preceding 
assumptions are : 

li+^es = t ies jUb + r ies E i+l,€s + pies 
pies = r les pies + ^ies jji+1, es + pies 

H 1 - h leo R 1+1 £a ) 


( 1 ) 

( 2 ) 

(3) 


L^)]o 

Mo 
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ales 


3 ies& 


(4a) 


vhere the range of 8 is from 0 to q and q is the number of 
generations of daughters. 


dpi®s0 ies0 icq -i+l.ea* -iesO _ies0 
— ^ — = nea(I +E ' ) - X P 


and for d = 1, 

^piesd 

dT 


K ies,d-1 piE£,d-l 


^iesd piesd 


(4b) 


(4c) 


The derivation of equations (4h) and (4c) is presented in the 
appendix. 

In the preceding equations : t^f s is the power-transmission 

coefficient for radiation of type s and energy € incident upon 
the ik* 1 station, which is transmitted with energy e; r^f 8 is 
the -power hack-scattering coefficient for radiation of type s 
and energy € incident upon the i^ station that is reflected as 
radiation- of energy e; h^ 6 ^ is the rate of conversion of energy 
of nuclear radiation of type a and. energy £ in the i th sta- 
tion of the absorber; is the rate of conversion of radi- 

ation of type a and energy £ absorbed in the i 01 station to 
radioactivity of a parent atom emitting radiation of type s and 
energy e; A* 68 ® is the decay constant for the radioactive parent 
isotopes in the i^ 1 station that emits radiation of type s 
and energy e; K^c^lrl the rate of conversion 'of radioactive 

energy in the i t:tl station of the (d-l)^ 11 generation of the iso- 
tope emitting radiation- of type Z and energy E to radioactive 
energy of the d^ generation of this isotope that emits radiation 
of type s and energy e. Capital Greek subscripts are used here 
to emphasize that only a single type and energy is indicated 

issd. 

because the increase in P is due entirely to the decay of 

one particular isotope. 


The interactions of radiations with the absorber that result 
in change of type are all contained in the radioactivity term P. 
The inclusion of all the interaction terms into P is physically 
plausible because a change of radiation type can occur only by 
means of a nuclear reaction. 
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In order to solve the problem, the set of equations (1) and 
(2) are need to express all the I's and E's in terms of the P's, 
and the boundary conditions [i-^T^o and ]jR n+ l( t)J g* The 
resulting expressions are then substituted in equation (4b) and 
the set of equations (4b) and (4c) becomes expressible in the fol- 
lowing form: 


dP 1 


es 


dT ’ 


A ies pies + pies [il( T )] 0 + C les [p n+1 (T)] 0 + D ies 


(5) 


If the indices are emitted, the solution to equation (5) may 
be written in the matrix form 

[p] = eL A l T |[p ] 0 + e- [A] T [f(t)] d T | (6) 


where 


[f(t)J = £b] [iHtiJq + [c] [B n+1 (T)j 0 + [i>j 


fA],[B],[0],[D] 


matrices all the elements of which are constants 


After all the P's have been found, I 1 , B 1 , and may be 

found as in reference 1. 


Example , - Assume that; 

(1) The absorber is three stations in thickness, 

(2) The initial incident radiation is 


: 1 1 7 


T 1 2 <x| 


T = 0 


T = 0 


(10-5 cos jtt) roentgen per hour' 


= (100-50 cos kt) roentgen per hour 


(3) The alpha-radiation induces a 1 Mev gamma- emitting iso- 
tope in the second station; this isotope decays to a 3 Mev neutron 
emitter which, in turn, decays to a stable isotope. 



8 


MCA W 2108 


(4) The 3 Mev neutrons induce 2 Mev gamma -rays in the source of 
incident radiation. 

(5) The other houndary conditions are as follows. At T » 0: ■ 

E 4 e t = H i t ji 2 y = p i o 0 

(6) The pertinent physical constants are; 

t 1 1 7 = t 1 2 7 = t 2 1 7 = 0.9 

t 2 2 7 = t 3 1 7 = 0.8 

t 3 2 7 = 0.9 

t 1 2 a = 0.8 

t 2 2 a = o.7 

t 3 2 a = 0.5 

t 1 3 n = 0.9 

t 2 3 n «= 0.8 

t 3 3 n = 0.9 

r 1 e y * 0 

r 1 2 a = 0.1 

r 2 2 a e 0.2 

r 3 2 a = 0.3 

r 1 3 n = 0.1 
r 2 3 n _ 0>2 

r 3 3 n = 0.1 
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K 2 3 n 1 = 0.5 second." 1 

■vO 2 a 0 ,0 2 7 0 , .-1 

A = A ' = 1 second 

A.2 1 7 0 = 10" 1 second" 1 
X 2 3 D 1 . 10 second”'*" 

H let = (i 3et = A let5 = 0 
^ 2 a = 10 ~ 3 second" 1 

3 n “ 10”^ second" 1 
All h = 0 


Solution . - The preceding data and conditions assume no 
Induced alpha -act ivity . The alpha-radiatlohs may therefore he found 
simply by solving the simultaneous algebraic equations involved. 

The methods outlined in reference 1 are suitable and when applied 
yield: 


jl 2 a 

JL 


10 1 



8 ‘ 8 

E 1 2 a 

— 

1 _8 1 _ 



co 

0 

H 

00 

_l 


The final results are: 


10 

7 




3 


T 4 2 a 

5 



5.9 


n 

10 50 


u 


I 2 2 a = 0.8296 I 1 2 a 

j3 2 a _ 0.6178 I 1 2 a 

I 4 2 a = 0.3089 I 1 2 a 
R 1 2 a = 0.3364 I 1 2 a 
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E 2 2 a = 0.2956 I 1 2 a 


E 3 2 a = 0.1853 I 1 2 a 


The neutron intensities may all "be expressed in' terms of the 
3 Mev activity of the second -generation daughter product in the 
second station. Thus, 

.j.4 3 n _ p2 3 n 1 
j3 3 n _ 10 p2 3 ‘n 1 

j23n_lp23nl 

9 

E 1 3 n p2 3 n 1 


E 


2 


3 n c 10 p 2 3 n 1 
9 


E 3 3 n 


1 p 2 3 n 1 
9 


The calculations of the intensities of the gamma-rays are more 
complicated; the methods of Case II (reference l)‘are applicable to* 
these computations and yield" the following results : 

I 2 1 7 = 0.9000 I 1 1 y 

I 2 2 7 = 0.9000 I 1 2 7 


: 3 1 7 a p2 1 7+ 0.81 I 1 1 7 

1 3 2 y = 0.7200 i 1 2 7 

1 4 1 7 = 0.8000 P 2 1 7 + 0.6480 I 1 1 7 
I 4 2 7 - 0.6480 I 1 2 7 

E i 2 7 = E 31 7 I= o 

E 1 1 7 = 0.9000 P 2 1 7 . 

E 2 1 7 . P 2 1 7 . 
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If P's are now determined, the problem will be solved. 

p2l7_p2l70 

The differential equation satisfied by p2 1 7 0 p s 
dP 2 ^ 1 7 0 = 10 -3 (j 2 2 a. E 3 2 a) . 1 p 2 1 7 0 

= HT 3 X 1.0149 il 2 a 1 p 2 1 7 0 


This equation is integrated and the given boundary conditions 
are substituted to result in 

+ it sin ict\ 

Similarly, 

n p2 3 n 1 

— a r= lp 2 l 70 _ lop 23 nl 

dT 2 

If the integrated expression for p 2 1 7 0 is substituted 
into this equation, it is found that 

p 2 3 n 3. = 0.05075 + 0.0002073 cos *T - 0.0007269 sin *T 
Finally, 

0 2 t 0 

dP ' = 10 -2 E 1 3 n _ p 0 2 7 0 

dT 

_ po -2 p2 3 n 1 m pO 2 7 0 


p 2 i 7 0 = p 0149 _ 0.005136 

V 10 


Therefore, 

p° 2 7 0 = 0.0005075 + 0.000002292 cos «T - 0.0000000696 sin jtT 

All the P's are to be evaluated at T = 1 day = 86,400 seconds. 
The values resulting therefrom are substituted into the preceding 
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expressions for the I*s and R's. Tn a n m u oh as 

jl 2 y _ p0 2 7 p0 2 y 0 , -fcR e following results are readily 

obtained : 

j2 2 a ^ 41,48 roentgens per hour 
I 3 2 a _ 30.89 roentgens per hour 
j4 2 a _ 15.45 roentgens per hour 
Rl 2 a = 16.82 roentgens per hour 
R 2 2 a = 14.78 roentgens per hour 
R 3 2 a _ 9,265 roentgens per hour 
j2 3 n _ 5.662 x 10" 3 roentgen per hour 

1 3 3 n c 5.662 x 10 ~ 2 roentgen per hour 

1 4 3 n = 5.0957 X 10 ~ 2 roentgen per hour 
r!3 n c 5.0957 x 10“ 2 roentgen per hour 
r 2 3 n _ 5,072 x 10“ 2 roentgen per hour 
r 3 3n _ 5.072 x 10 “ 3 roentgen per hour 
j2 1 7 _ 4.500 roentgens per hour 

1 2 2 7 = 4.588 X 10 “* roentgen per hour 

1 3 1 7 = 5.064 roentgens per hour 

I 3 2 7 = 3.670 x 10 “ 4 roentgen per hour 

1 7 = 4,052 roentgens per hour 

2 7 _ 3.303 x 10" 4 roentgen per hour 
R 1 1 7 = 9.130 X 10 -1 roentgen per hour 
r! 2 7 r= 0 roentgen per hour 

R 2 •*- y = 1.015 roentgens per hour 



MCA SET 2108 


13 


g2 2 7 = o roentgen per hour 
.R 3 1 7 = 0 roentgen per hour 
r 3 2 7 - o roentgen per hour 


DISCUSSION 

Many induced radioactiTe isotopes decay to one or more inter- 
mediate radioactivities "before becoming stable atoms. Because, in 
many cases, the radiations from these intermediate isotopes are 
more harmful than the original activity, the entire gamut of radi- 
ations must be accounted for. 

In addition, the possibility exists that radioactivity from 
the absorber will induce a radioactivity -in the source itself that 
differs from the original radiation, theoretically, this consider 
ation may be as important as the one in the preceding paragraph 
because the decay of almost all the induced radioactivities can 
be predicted from the available data, whereas the effect of the 
back-reflected radiation on the source is not as readily obtained. 
Furthermore, any change in the source results in possible subse- 
quent changes in the radiation at the various stations of the 
absorber. 

The methods of references 1 and 2 have been generalized to 
include the preceding two considerations. The main limitations of 
the method are: 

« 

(1) The assumption of normal incidence at all stations of the 
absorber 

(2) The assumption of stations of the absorber as plane 
parallel surfaces 

(3) For any practical problem, a large number of physical 
constants are required for solution. A great many of the 
required constants are not yet available in the literature. 
Nevertheless, the method should currently prove useful for 
obtaining approximate results in absorption problems in 
which reasonable estimates of the constants may be made. 
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SUMMARY OF RESULTS 

Tlie problem of determining transmission and absorption of time- 
dependent radiation through thick absorbers has now been generalized 
to include the following time -dependent cases: 

(a) Radiation is normal to an absorber of which the stations 
are plane parallel surfaces. 

(b) An arbitrary number of energies and types of radiation 
are involved in the source and in the radioactivities of the 
absorber. 

(c) Radioactive isotopes decay to stable atoms in a finite 
number of steps . 

(d) The source is an arbitrary function of time and its activ- 
ity may be affected by the back-reflected radiation from the 
absorber. 

Thus, the amount of heat and the activities in various portions 
of an absorber can be calculated for problems of complexity embody- 
ing cases (a), (b), (c), and (d) . 


Lewis Flight Propulsion Laboratory, 

Rational Advisory Committee for Aeronautics, * 
Cleveland, Ohio, December 9, 1949. 
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APPENDIX - DERIVATION OP EQUATIONS 4(b) AND 4(c) 


Let N a iesd equal the number of atoms in the station 

of the absorber in the generation that are capable of emit- 

ting radiation of type s and energy e and let the number of 
particles of radiation be denoted by jfi® 8 ^. Ehe following equa- 
tion may then be written: 


d jjiesd _ Qiesd _d_ ^ iesd 

dT dT a 


(Al) 


where C iesd is a constant that accounts for the possibility that 
in the decay process an atom may emit more than one radioactive 
particle of type s and energy e. 


The following two equations define the rate of change of 
jt iesd. 

"a 


a ies 

± N fl les0 = + b 1+1 > €CT ) -A ies0 N ies0 

dT Ei SRo 


(A2) 


where the 0 generation is defined as the original induced member 
of the radioactive decay series. The term a 1 !^ is the fraction 
of the energy of radiation of type O and energy e absorbed in form- 
ing radioactive atoms of N a ies0 and E^sEq i0 ener &y required 
to change a stable atom to the radioactive parent of the series. 

For d=l 

N a iesd = A ies > d “ 1 N a ies * d ” 1 - X iesd N a iesd (A3) 

dr a a a 

Now 


pies _ pies 5 
P P 5 

where the range of 5 is 0 to q. 
However , 


piesd = jji 




Qiesd ^iesd jj iesd 
& 


(A4) 


(A5) 
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where E^p p ia the energy emitted, ■when a radioactive atom of 
d d+1 

the d'kk generation changes to a radioactive atom of the 
(d+l)^* 1 generation. Therefore, 


dP 


lead 

dT 


E- 


K,R. _ 
d d+1 


c iead ^iesd 




iesd 


(A6) 


If equations (A2) and (A3) are substituted into equation (A6), 
there results : 


_ pi c iesO j. iesO 
dT EoEl 0 x 


a 1 " (I l€ ° . B 1+1 ’ €a ) 


E 


SR r 


^iesO N iesO 


(A7) 


and for d e 1, 
3 iesd 


dP- Lt ’ ua _ „iesd -v iesd 


dT R,B, , 

d d-1 


C X ■ 


^ies, d-1 


1 ies, d-1 _ 


a 


^iesd jj^iesdj 


(AS) 


Finally, if equation (A5) is substituted into equations (A7) 
and (AS), they may be written 


dP 


,iesO 


_ ^iesO (jiea + * R i+l,ea) _ ^iesO piesO 


(as) 


and for d = ly 


d-P iead = K ies, d-1 piES, d-1 _ ^iesd piesd 


dT ES 


(A10) 


where n^- ea ® is 
€G 


i?i niesO -viesO „ies 
E BqRp 0 A a eo 


E' 


SRr 


o 
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and K 3 -®^ (i “ 1 is 
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